Onset of an Insulating Zero-Plateau Quantum Hall State in Graphene 
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We analyze the dissipative conductance of the zero-plateau quantum Hall state appearing in 
undoped graphene in strong magnetic fields. Charge transport in this state is assumed to be carried 
by a magnetic domain wall, which forms by hybridization of two counter-propagating edge states 
of opposing spin due to interactions. The resulting non-chiral edge mode is a Luttinger liquid of 
parameter K, which enters a gapped, perfectly conducting state below a critical value K c ~ 1/2. 
Backscattering in this system involves spin flip, so that interaction with localized magnetic moments 
generates a finite resistivity R xx via a "chiral Kondo effect". At finite temperatures T, R X x(T) 
exhibits a crossover from metallic to insulating behavior as K is tuned across a threshold K mi ■ For 
T — > 0, Rxx in the intermediate regime Kmi < K < K c is finite, but diverges as K approaches K c . 
This model provides a natural interpretation of recent experiments. 

PACS numbers: 71.10.Pm, 72.10.Fk, 73.20.-r, 73.23.-b, 73.43.Nq 



Introduction- Graphene, a honeycomb network of 
carbon atoms, is perhaps the most remarkable two- 
dimensional system to be studied in the last few years [l[ . 
This system differs from the standard two dimensional 
electron gas (2DEG) in supporting two Dirac points in 
its fermion spectrum which are at the Fermi energy when 
the system is nominally undoped. In the presence of a 
magnetic field, the Landau level spectrum differs from 
the standard 2DEG in supporting positive and negative 
energy (Landau level) states, as well as a zero energy 
(lowest) Landau level in each of its two valleys. In the 
earliest experiments, this last property was understood 
to account for the absence of a quantized Hall effect at 
filling factor v = 

In stronger fields, and with higher quality samples, a 
plateau does appear to emerge near v = which is likely 
associated with a resolution of the spin-splitting of the 
lowest Landau level 0, H, [H] . However, the behavior of 
the longitudinal (dissipative) resistance (Rxx) at this fill- 
ing has been difficult to explain: different measurements 
have demonstrated that it may either decrease [B} or in- 
crease @ with falling temperature T in samples that ap- 
pear to be quite similar p}. Moreover, in Ref. Q the 
v = quantized Hall state exhibits a transition to an in- 
sulator at a critical magnetic field H c , where R X x(T — > 0) 
diverges. The scaling of R xx as a function of field H be- 
low H c appears to signify a quantum phase transition of 
the Kosterlitz-Thouless (KT)@| type. 

In this paper, we demonstrate that such a variety of 
transport behaviors can be explained quite naturally if 
the system contains local magnetic moments near its 
edge, which may occur due to chemical passivation of 
dangling bonds or imperfections in the lattice Inter- 
action with such magnetic degrees of freedom provides 
a mechanism for back-scattering in the primary channel 
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FIG. 1: (Color online). Phase diagram of the DW coupled to 
a magnetic impurity via a fixed Kondo interaction. The solid 
blue curve depicts the gap A 3 (K). The dashed lines mark the 
boundaries between distinct transport regimes (see text). 



for charge conduction at v = 0, a magnetic domain wall 
(DW) EH EH which may carry current that is unaf- 
fected by ordinary, non-magnetic impurities. Transport 
properties of this mode are dramatically altered by tun- 
ing of a Luttinger parameter K [l3| , which in this system 
depends on the details of the edge potential and on the 
magnetic field H. 

Our main findings are summarized in the phase dia- 
gram depicted in Fig. [TJ For K > K c , the ideal system is 
a Luttinger liquid, with K large enough that a Kondo im- 
purity renders the system insulating at zero temperature 
[Rxx — * oo]- As K is tuned below K c , the impurity- free 
DW system undergoes a quantum phase transition to a 
perfect conducting state protected by a gap A s (K) (solid 
line in Fig. 1). In this case a Kondo impurity generates 
a finite resistance at zero temperature, approaching this 
value from below (dR xx /dT < 0). The unusual thermal 
behavior of this "odd metal," results from a competi- 
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tion between two relevant perturbations which separately 
would render the system either perfectly conducting or 
truly insulating. Finally, for K < Kmi, a Kondo impu- 
rity is an irrelevant perturbation, yielding more typical 
metallic behavior, dR xx /dT > 0. This variety of behav- 
ior is consistent with existing published data [B|, @ • 

The remarkable edge-mode structure of the DW can 
already be seen within a non-interacting picture. The 
lowest Landau levels support an electron-like and a hole- 
like mode for each spin, which cross precisely at the Fermi 
level when Zeeman splitting is introduced [see Fig. [5] 
Tot 11 , 14 1 . Accounting for interactions within a Hartree- 
Fock description, the counter-propagating states admix 
to form the DW [13], with an unusual collective mode. 
In particular, due to the spin-charge coupling inherent 
in a state projected into a single Landau level [HI, it 
supports topological charged excitations which can prop- 
agate freely along the edge. In addition, a spin gap may 
open in the collective excitation spectrum when interac- 
tions are not too strong, with important consequences for 
dissipative transport. 

A prominent candidate for a dissipation mecha- 
nism in this system is the scattering between counter- 
propagating edge states [Tf|, which one might expect 
to be enhanced by the DW structure which strongly 
hybridizes forward and backward propagating states. 
However, the requirement of spin conservation forbids 
backscattcring within a single DW by static impurities 
111 ]. In contrast, scattering from local magnetic impuri- 
ties does allow backscattering and generates dissipation. 
As we show below, the coupling of a DW to a magnetic 
impurity leads to an unusual "chiral" Kondo effect, since 
left and right moving electrons each come in only a single 
spin flavor. The resulting resistance at finite T exhibits 
the variety of behaviors summarized above, which we dis- 
cuss below in more detail. 

Model - As a first stage we derive an effective model 
for a clean graphene sheet, with a straight edge along 
the y direction. Our starting point for the analysis is 
to view the system as a two-dimensional (2D) quantum 
Hall ferromagnet, in which only the two Landau levels 
closest to zero energy are included. These levels carry 
opposite spin, and may be grouped together into a spinor 
whose long- wavelength degrees of freedom are those of a 
Heisenberg ferromagnet, with Hamiltonian T^, HI 



H Q 2D = I d 2 r- 



J 



£|9 M S(r)| 2 + A(x)S z (r) 



(1) 



Here J represents an exchange stiffness due to interac- 
tions. A(x) encodes the combined effects of the Zee- 
man interaction and the edge potential (dictated by 
the boundary conditions and electrostatic environment 
fl7j]): it increases monotonically from the bulk (neg- 
ative) Zeeman energy for x well inside the bulk to a 
positive value at the edge, leading to the energy cross- 




FIG. 2: (Color online). Non-interacting energy levels near 
the edge. Solid black lines are spin up states, dashed red lines 
are spin down states. The dotted line denotes the chemical 
potential at zero energy. 



ing in the non-interacting spectrum illustrated in Fig. 
[2] Initially treating the system classically, we parame- 
terize the spin field by S(r) = Sil(r), with S = 1/2 
and Q = fsin^(r) cos y(r), sin^r) sin y(r), cos^(r)]. As 
shown in [121 ]. the form of A(x) in the above energy func- 
tional dictates a minimum energy configuration where x 
is spatially constant, and tp(r) = ^o( x ) exhibits a non- 
trivial topology: tpo(x — ► — oo) — * 0, ipo(x — ► edge) — > n. 
A domain wall (DW) is thus formed parallel to the edge. 
The ground state energy is independent of x, implying 
a broken symmetry and an associated collective mode 
propagating with a wavevector k y along the DW [lU . 

The 2D semiclassical theory may be projected into 
its low-energy subspace, containing the one-dimensional 
(ID) mode. This theory is identical in form to a semi- 
classical treatment of an XXZ antiferromagnetic spin-1/2 
chain; an appropriate choice of coefficients for the latter 
allows the models to coincide at the semiclassical level. 
As the spin chain incorporates the spin-1/2 nature of the 
system, we work with this model, with lattice spacing set 
by the magnetic length t = ^Jhc/eH : 



Hd\V — Jxy ^ [SjSj + 1 



(2) 

Here J xy is related to the exchange energy J via J xy — 
JS 2 Af, with M = J2 X sm2 ^o(x) (ipo(x) the classical 
DW texture), and J z depends on the edge potential: 
Jz = jpzs T,x A(x)sin 2 V>o(20[l ~ cos-0o(x)]. Using a 
standard Bosonization scheme for the spin operators in 
the continuum limit II 311 - 



S±(y) 



V2t 



[(-) y + cos(20)] 



dy<j> , (-)* 



7T 7ra 



cos(2</>) 



(3) 



we map How to a (single-flavored) Luttinger model with 
a sine-Gordon correction. In units where H = 1, 

H DW = I ^ \uK(irtt) 2 + ^-(dyef)) 2 ~ ^ cosm)} 
J Ztt y K ira z J 

(4) 
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where the Boson field <j)(y) and IL(y) = are canoni- 
cally conjugate. The spin- wave velocity u and Luttinger 
parameter K are related to the XXZ parameters by 



uK Jxy i j^- — Jxy 



IT 



(5) 



and the short distance cutoff a ~ I. 

The effective model for the ideal DW system [Eq. (QJ] 
exhibits a quantum phase transition upon tuning of the 
parameter K across a critical point K c , from a Luttinger 
liquid (LL) phase at K > K c [where the cos(4</>) term 
is irrelevant] to an ordered [gapped] phase below K c . In 
terms of the spin fields, the latter phase is characterized 
by ordering of the field S z , and the opening of a gap 
A s to spin-flip excitations. The critical behavior near 
K c can be derived from a perturbative renormalization 
group (RG) analysis [3], yielding a ifT-transition at 
K c ~ 1/2. As K approaches K c from below, the gap 
tends to vanish according to the scaling law 



suppressed for T < A s . Its 'perfect conduction' is hence 
more robust. Indeed, Hdw [Eq. ©] is analogous to a 
low-energy model of a quasi ID superconducting system, 
where two parallel wires are Josephson coupled. The field 
4> describes the relative superconducting phase. u/K the 
superfluid stiffness, uK the charging energy [19J, and the 
last term in How is the inter-wire Josephson coupling. 
There as well, a local perturbation analogous to Hk [Eq. 
(0)] is required to generate finite dissipation via a phase- 
slip mechanism. This occurs when a Josephson vortex 
(i.e., out-of-phase current configuration) penetrates into 
the bulk of the double- wire system. In the gapped phase, 
where Josephson coupling is relevant, such processes are 
suppressed beyond the Josephson length (~ 1/A S ). 

We next consider a finite but small Kondo interaction 
and T > A s , and calculate the dissipative part of the 
electric conductance to leading order in Hk' 



A S (K) 



■ exp 



C 



-1 f L/2 

<5G=lim— / dy^{x(y,y';u)} 

u^-Q LUJ J_ L / 2 



(10) 



(6) 



with C a constant of order unity [see Fig. [T]. 

We now assume that a local moment at the origin (rep- 
resented by a spin- 1/2 operator (?) couples to the original 
2D spin system via an exchange interaction JS(r — 0) •<?. 
In order to project this into the low-energy subspace, we 
express it in a rotated basis where the local moment is 
coupled to fluctuations about the ground state DW con- 
figuration. This yields an effective Kondo coupling to the 
ID spin chain operators: 



where \ — (je(y)',je{y'))Z is the retarded correlation 
function evaluated with respect to Hq (see, e.g., fl3j| ) . 
and L the length of the system, assumed finite. Note that 
SG is a negative correction to the ideal conductance, so 
that \SG\ is proportional to the backscattering rate and 
consequently to the longitudinal resistance R xx [13] • Us- 
ing the bosonized representation Eqs. [[3"]).(j3]) and the 
expression for the local spin correlator (for t > 0) 
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(ii) 



Hk = -|p-(5+(0)cr_ + S-(0)a + ) + J K S z (0)a z . (7) we obtain the T-dependent conductance 



The resulting Hamiltonian for the effective ID system 
then takes the form H = Hq + Hk , with 



H = H 



DW 



(8) 



h \ u 



Here gx oc J K y is dimensionless, k = 1/2K and 




(12) 



Note that E z is an effective Zeeman energy of the local- 
ized spin, associated with a local magnetic field which 
depends on its position within the width of the DW. 

Transport - To derive the transport properties of the 
system, we first express the charge current operator in 
terms of the DW fields. Using the standard definition 

S Hpw 



G(z,e) 



sin 



(?) 



(fli- 



rt I \ f nK \ 

2m[z) cos y~2~ ) 



t cos(zt) 
sinh(7rf)] K 

t sin(^i) 



dt- 



sinh(7rf)] K 



(13) 



J 



SA , with A(y,t) a vector potential, we find 



3e{y,t) = 2eJ z S z (y,t) 



(9) 



This remarkable relation between spin and charge cur- 
rent encodes a unique property of the DW: in the ab- 
sence of spin- flip processes (i.e. for J K V — 0), current is 
conserved and the DW behaves as a perfect conducting 
channel. The distinction between its two phases becomes 
apparent when spin impurities are added: in the gapped 
phase, spin-flip and hence back-scattering excitations are 



with m(z) = [f(z) — /(— z)] the local spin magnetization. 

In the finite T regime where T ^> E z , the leading In- 
dependence of SG(T) inferred from Eqs. (fT2|) . (fT3|) is a 
power-law, which changes sign at k = 2: 



5G{T) 



(14) 



The conductance therefore exhibits a crossover from a 
metallic behavior (dG/dT < 0) for n > 2 to an insulating 
behavior (dG/dT > 0) for k < 2. This crossover would 
be a true metal-insulator (MI) as K is tuned through 
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Km i ~ 1/4, as depicted in Fig. ([TJ) , if not for the exis- 
tence of the gap. Note that K mi is perturbatively renor- 
malized to values below 1/4 by and J^-, so that one 
always finds Kmi < K c . The MI crossover then always 
occurs within the gapped phase. Thus the power-law 
behavior of G{T) is restricted to finite T ^> A s . 

We now focus on the intermediate regime Kmi < K < 
K c . Assuming further that E z ~ 0, we consider the low 
T limit T <C A s where the field cf> is ordered, while 9 is 
strongly fluctuating. The power-law decay of the corre- 
lation function (^e l9 ^e~ ie ^) (characteristic of the LL) 
is modified by a rapidly oscillating exponential factor: 

<■""'«- w, >"~©M-1r + KK)} 

(15) 

where the decoherence rate 1 = A 2 s L/2-ku corresponds 
to a 'condensation energy' in the superconductor ana- 
logue. The resulting 8G saturates as T — > and becomes 

SG ~ -t\~ k . (16) 

It then follows from the definition of t$ and Eq. ([6]) that 
the resistance diverges when K approaches K c as 

"■--"-'{tpi}- (17) 

Here C = (7(2 - k); for k - l/2Jf c ~ 1, C ~ C. 

We finally comment on the possible relation of our re- 
sults to the experimentally observed R xx in the zero- 
plateau state. Provided K is tunable by a magnetic field 
H, the system is shown to exhibit diverse transport prop- 
erties, which could explain the apparent discrepancy be- 
tween the data of Refs. [5J] and [6j. Furthermore, we find 
a quantum critical point associated with the closing of a 
spin-gap, manifested by a divergence of R XX (T — > 0) ac- 
cording to a scaling law characteristic of a KT-transition 
(in 1 + 1 -dimensions). The latter is remarkably reminis- 
cent of the peculiar field-tuned transition to an insulator 
observed in [|| . A rough estimate of K in terms of phys- 
ical parameters yields K ~ \J e eA > w ith A av an aver- 
age edge potential strongly dependent on sample details. 
In particular, when A av is dominated by electrostatic ef- 
fects induced by the gate voltage [13] , it is typically larger 

2 

than the exchange energy |j. This would yield K < 1 
and monotonically increasing with H , as assumed by our 
theory. 

Another assumption of our theory which requires jus- 
tification is the neglect of E z . A finite E z introduces 
a cutoff which competes with the gap A s , and a finite 
magnetization of the local moment which reduces its con- 
tribution to backscattering. We note, however, that the 
realistic system presumably contains a multitude of spin 
impurities with randomly distributed local parameters. 
The dissipative resistance is then dominated by the mo- 
ments with minimal E x , and our analysis applies as long 
as it is smaller than all other energy scales. 



To summarize, we study a model for charge transport 
carried by fluctuations of a DW propagating along the 
edge of a graphene sample in the v — quantized Hall 
state. Backscattering induced by localized magnetic im- 
purities at the edge provides a dissipation mechanism. 
Its competition with a spin-gap generated in the ideal 
DW may possibly explain the rich variety of conductance 
properties observed in recent experiments. 
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